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Abstract 

We summarize the science opportunity, design elements, current and projected partner ob- 
servatories, and anticipated science returns of the Astrophysical Multimessenger Observatory 
Network (AMON). AMON will link multiple current and future high-energy, multimessenger, 
and follow-up observatories together into a single network, enabling near real-time coincidence 
searches for multimessenger astrophysical transients and their electromagnetic counterparts. 
Candidate and high-confidence multimessenger transient events will be identified, character- 
ized, and distributed as AMON alerts within the network and to interested external observers, 
leading to follow-up observations across the electromagnetic spectrum. In this way, AMON 
aims to evoke the discovery of multimessenger transients from within observatory subthresh- 
old data streams and facilitate the exploitation of these transients for purposes of astronomy 
and fundamental physics. As a central hub of global multimessenger science, AMON will also 
enable cross-collaboration analyses of archival datasets in search of rare or exotic astrophysical 
phenomena. 
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1. Introduction 

Wc stand at the dawn of multimessenger astrophysics - a quest to use the messenger 
particles of all four of nature's fundamental forces to explore the most violent phenomena in the 
universe. Observatories first imagined a generation ago are finally being realized, including the 
Advanced LIGO [1] and Virgo [2] gravitational- wave detectors, the ANTARES [3] and IceCube 
[4] high-energy neutrino observatories, and the Pierre Auger Cosmic Ray Observatory [5]. On 
the ground and in space, they are complemented by high-energy observatories including the 
Swift [6] and Fermi [7] satellites, the HESS [8], VERITAS [9], and MAGIC [10] TeV gamma-ray 
telescopes, and the HAWC [11] TeV gamma-ray observatory. 

Collectively, these facilities promise the first detections of gravitational waves and high- 
energy cosmogenic neutrinos, the resolution of the mystery surrounding the origins of ultrahigh- 
energy cosmic rays, and a new window into the formation and evolution of black holes. Given 
their nature as first-generation facilities, however, the sensitivities of the non-electromagnetic 
observatories are naturally limited, with rates of detection for transient events of publishable 
significance known or expected to be low, perhaps a handful per year (or, in the case of 
Advanced LIGO and Virgo at design sensitivity, a few dozen [1]). 

During the intervals prior to and between detection of these rare high-significance events, 
the multimessenger facilities will be buffeted by signals from a far greater number of lower- 
significance events that will be statistically indistinguishable from background or noise pro- 
cesses. Such "subthreshold" events are, by definition, unrecoverable as astrophysical signals in 
the data stream of any individual faciUty. However, if they arc accompanied by a subthreshold 
signal in another multimessenger channel they can be identified, and potentially achieve high 
significance, via careful coincidence analysis of the data streams from multiple facilities. 

In this paper we present the scientific case for the Astrophysical Multimessenger Obser- 
vatory Network and describe its important elements. AMON will weave together existing 
and forthcoming high-energy astrophysical observatories into a single virtual system, capa- 
ble of sifting through the various data streams in near real-time, identifying candidate and 
high-significance multimessenger transient events, and providing alerts to interested observers. 

As we show, AMON will enable a significant enhancement in the effective aggregate sensi- 
tivity of the world's leading multimessenger facilities for a small fraction of the facilities' total 
cost, provide the first near real-time alerts for multimessenger transient sources, and simplify 
the mechanics and politics of cross-collaboration analyses for all partners. As such, we believe 
AMON represents a natural next step in the extension of the global astronomical community's 
vision beyond the electromagnetic (EM) spectrum. 

The development of AMON is currently underway. Signatories to the AMON Memoran- 
dum of Understanding^ (MOU) include the IceCube and ANTARES neutrino observatories, 
the HAWC TeV 7-ray observatory, and the Swift orbital telescope. Exploratory discussions 
with the Pierre Auger Cosmic Ray Observatory, VERITAS, LIGO (including the GEO-600 
collaboration), and Fermi scientific collaborations have led to signed letters of commitment, 
with negotiations toward MOU signatures from all parties ongoing. Discussions have also been 
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initiated with candidate follow-up facilities including ROTSE-III [12] and the Palomar Tran- 
sient Factory [13], with the goal of bringing these observatories into the collaboration prior to 
or shortly after the commencement of real-time AMON alert operations. 

AMON is structured as an open and extensible network, with an MOU that allows straight- 
forward incorporation of new triggering and follow-up facilities. Collaborations interested in 
the scientific goals of AMON, and with useful triggering or follow-up capabilities to contribute, 
are encouraged to contact the authors for information about joining AMON. First versions of 
the necessary supporting hardware and software infrastructure for AMON are being installed 
at Penn State, and initial analyses on archival and simulated real-time data streams will get 
under way shortly thereafter, as a means of preparing to bring the first set of triggering facilities 
on-line within the next year. 

The paper is structured as follows: In Section 2 we provide an overview of the scientific 
opportunity for AMON which motivates our efforts. In Section 3 we discuss the elements 
of AMON, including the technical and operational protocols that we propose to adopt in 
linking the partner facilities, and the algorithms that we will use to identify coincident signals. 
Section 4 presents detailed simulations of multimessenger transient sources, and additional 
theoretical case studies, which illustrate the gains that stand to be realized by AMON. Section 5 
presents our summary and conclusions. 

2. The AMON Science Opportunity 

AMON is intended to contribute in several ways to the first decade of multimessenger 
astronomy. A common focus of these approaches is on multimessenger transient events that 
are observed as coincident (potentially subthreshold) signals in the data of AMON partner 
facilities corresponding to two or more distinct types of messenger particle. In this section, 
wc briefly review the strongest candidates for these multimessenger transient sources and the 
current theoretical expectations for their properties, rates, and broader implications for physics 
and astrophysics. 

2.1. High-Luminosity Gamma-Ray Bursts 

As the most violent cataclysms known, and as the sites of the most highly-relativistic 
cosmic outflows, gamma-ray bursts (GRBs) have long been considered likely multimessenger 
transient sources. Recent progress in distinguishing the distances and energetics of GRBs and 
understanding the likely nature of their progenitors [14] means that we can review the expected 
multimessenger signals of the different varieties of GRB on a case by case basis. 

The traditional or "classical" long-duration, high-luminosity gamma-ray bursts (HL-GRBs) 
are the most common type of GRB detected by satellite experiments, being observed as bright 
seconds- to minutes-long bursts of 7-radiation from high-redshift, z ^1. HL-GRBs arc believed 
to arise when a massive star (M ^ 25Mq) undergoes core collapse to a black hole (BH); 
confirmation of this "collapsar" model [15] for the HL-GRBs has been most dramatically 
provided by spectroscopic observations of ensuing type Ibc supernovae (SNe), with other lines 
of evidence also contributing [16]. 

In the collapsar model, formation of a high angular-momentum BH and accretion of resid- 
ual gas through an accretion disk produces a relativistic jet. In "successful" GRBs the jet 
pierces the stellar envelope, accelerates to high Lorentz factor, and radiates gamma-rays for 
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tens of seconds, providing a bright electromagnetic trigger for observers within the jet colli- 
mation angle. The typical HL-GRB jet energy of ~ 10"''^ to 10^^ erg, collimated within an 
angle of {6) ~ 5°, roughly 1/250 of the sky, yields the observed isotropic-equivalent energies 
of Eo = E (47r/flj) ~ 10^^ to lO^'^ erg [17]. For a burst duration T ~ 10 s, this corresponds to 
an isotropic-equivalent luminosity L-y ~ Eq/T = 10^^ to 10^^ erg s~^. Given a 250:1 jet colli- 
mation factor, the nearest off-axis bursts are anticipated to be located 250^/^ ~ 6 times closer 
than the nearest on-axis bursts. Off-axis bursts, defined as bursts whose bright initial 7-ray 
emission does not illuminate observers at Earth, may still be observed via their less luminous 
shock-breakout emission, their prompt gravitational wave (GW) signal, or their subsequent 
afterglow and/or supernova components. 

In the standard internal-shock model, the gamma-ray emission of HL-GRBs is produced 
by internal shocks at a dissipation radius Vd F'^cdt ~ (3 x lO-*^^ cm) 77^ 5 (5i_2 [18], where 
?72.5 = r/300 is the bulk Lorentz factor divided by 300 and St-2 = dt/10~'^ s is the variability 
time of the central engine in hundredths of a second. These shocks Fermi-accelerate electrons, 
which produce gamma-rays via synchrotron and inverse Compton scattering, and are then 
boosted by the bulk motion of the relativistic outflow. For a discussion of alternative models 
see ref. [19]. 

The same shocks responsible for electron acceleration and 7-ray emission should also ac- 
celerate protons, leading to photo-produced pions which in turn yield high-energy cosmogenic 
neutrinos (CN) and 7-rays from charged and neutral pion decays, respectively [20, 21]. Initial 
assumptions were that the energy in relativistic protons would be of the same order as the 
energy emitted in 7-rays, Ep ~ [1/ f(,)Ej, where /g ^ 1 is the fraction of proton energy given to 
electrons (and observed as GRB photons), so that the optical depth Tp^ determines the CN and 
TeV gamma-ray luminosity, LtcV ~ Tp^Lp. The CN to 7-ray flux ratio expected in this model 
has been quantified [22] and used by IceCube [23-25] to set limits which are already a factor 
of five below the naive predictions. However, careful consideration of the underlying physics 
[26-28] suggests that IceCube observations will need to continue for several more years before 
the somewhat reduced CN fluxes of more realistic models are tested. Ultimately, detection of 
GRB-related CN is expected if HL-GRBs contribute a significant fraction of the highest-energy 
cosmic rays. 

Regardless of the total energy release of the HL-GRBs, which may well be an order of 
magnitude or more above the beaming-corrected 7-ray energies E ^ 10^^ erg, no significant 
GW emission will be produced if the core collapse, jet production, and burst processes maintain 
approximate axisymmetry throughout. However, it is possible for a rapidly rotating core 
and accretion disk to develop bar and/or fragmentation instabilities which could result in 
substantial GW emission. In the most optimistic case [29] the resulting GW signal is periodic 
and roughly as strong as the signal from neutron star binary mergers, and hence, visible out 
to similar distances (hundreds of Mpc for the advanced facilities) using ground-based GW 
detectors. 

While some 3D simulations of stellar collapse suggest a fragmentary process, and a cor- 
respondingly weaker GW energy release Eqw ~ 10~^Mqc^ ~ 10^^ erg [30], these have not 
included relativistic effects that are known to contribute to large bar asymmetries [31], as 
manifested in other simulations [32]. When present, these instabilities result in GW emission 
comparable to that found in optimistic analytical calculations [29]. 
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2.2. Low-Luminosity Gamma-Ray Bursts 

Low-luminosity GRBs (LL-GRBs) are underluminous long-duration GRBs, having substan- 
tially lower isotropic-equivalent energies, Eq ~ 10'^^ to 10^^ erg, than the HL-GRBs. Because 
of their lower 7-luminosities, they are currently only detected at low redshift, z ^ 0.5, where 
they provide the bulk of the observed GRB-SN Ibc associations [16]. Estimates of the LL-GRB 
rate suggest that they occur at ^100 times the rate of HL-GRBs while constituting ^1% the 
rate of type Ibc SNe overall [33-36] . 

The 7-ray emission of LL-GRBs may be due, in some or all cases, to a relativistic shock 
breakout [37, 38] and hence may not require the Lorentz factor F ^ 100 jet needed to explain 
the high-luminosity, high-variability 7-emission of HL-GRBs. Such shock breakout emission 
would also be expected to be uncollimated (uniform over the sky) or nearly so. 

Predictions for the CN emissions of LL-GRBs have been explored by [39, 40] for relativistic 
jet models and by [41] for shock breakout models. In relativistic jet models, a straightforward 
approach scales CN fluxes from HL-GRB predictions according to their 7-ray luminosities, and 
anticipates similar spectra with E'^(j)i, peaking at PeV to EeV energies. Given current evidence 
for less-relativistic jets in LL-GRBs, however, these predictions may be optimistic. 

Shock breakout models [41], by contrast, predict softer neutrino spectra peaking in the 
TeV to PcV range, with luminosities for reference events that would make individual bursts 
detectable to IceCube within a horizon of D ^ lOMpc. The LL-GRB rate within this horizon 
is thought to be small, r ~ 0.002 year~^, but prospects for detection could be improved by 
stacking analyses using ^dozens of more distant events identified by their prompt high-energy 
EM signature. 

Similarities in explosion energy, ejecta velocity, and synthesized nickel mass for the type 
Ibc supernovae of LL-GRB and HL-GRB events suggest that the details of core collapse for 
the two event classes, and their GW emissions, may be similar. However, given the critical 
role that high values of core angular momentum are thought to play in powering HL-GRB jets, 
and the need for accretion or disk instabilities to power significant GW emission, the relative 
weakness of LL-GRB jets may suggest correspondingly reduced prospects for GW emission. 
If the GW emissions are competitive with those of binary neutron star-neutron star (NS-NS) 
mergers in even some cases, this would make the LL-GRBs with their relatively isotropic 7-ray 
emissions a highly-promising target population for ground-based GW detectors. 

2.3. Short-Hard Gamma-Ray Bursts 

Short-hard gamma- ray bursts (SHBs), apart from their shorter durations and somewhat 

harder spectra, are observationally similar to HL-GRBs. Their harder spectra may indicate 
higher bulk Lorentz factors [18], while their reduced durations, tf, ^ 2 s, suggest shorter accre- 
tion timescales compared to those of the HL-GRB collapsars. A softer "extended emission" 
episode lasting for t ~ 30 to 100 s after the burst itself is present in about a third of SHBs, 
accounting for 5-20% of the total energy [42]. 

Thanks to afterglow discoveries and host galaxy identifications of the Swift era, consensus 
now holds that SHBs are likely due to compact binary (NS-NS or NS-BH) mergers [43, 14]. 
If so, these events are associated with progenitor systems that are observed within our own 
Galaxy (as relativistic NS-NS pulsar binaries), have merger rates that can be estimated from 
population synthesis modeling (e.g., [44]), and should emit strong GW emission in a highly- 
calculable "chirp" waveform across the frequency range of ground-based GW observatories 
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just prior to merger. Indeed, binary mergers are the primary extragalactic target population 
for next-generation LIGO and VIRGO, with expected event rates at fuh design sensitivity 
expected to reach dozens annually [1]. 

Afterglow observations for a handful of Swift-detected SHBs show evidence for typical 
collimation angles (0) ^ 6°, corresponding to beaming fractions of 200:1 (similar to that for 
HL-GRBs) or less [45]. These beaming corrections adjust the observed isotropic-equivalent 
energies of the bursts, Eq ^ 10^^ erg to lower inferred jet energies E ~ 10^^ to lO^'' erg. 
As such the nearest off-axis mergers are predicted to be observed at 200^^'^ ~ 6 times closer 
distance than the nearest on-axis events, where they may be detected as GW inspiral signals 
without SHE counterparts. 

Because of their very similar 7-ray characteristics and inferred jet properties, it is assumed 
that the radiation physics of the HL-GRBs and SHBs are similar, apart from possible slight 
differences due to the SHBs' reduced durations and potentially larger Lorentz factors and 
collimation angles [18]. CN production models may thiis conservatively be carried forward by 
assuming similar CN to 7-ray flux ratios as for the HL-GRBs. 

In terms of GW emission, SHBs in the compact binary merger model represent a "dream 
scenario" , as systems in which the GW waveform is known to be strong (detectable to hundreds 
of Mpc by advanced detectors) and calculable, and moreover, has already been implemented in 
data analysis systems via matched-filter algorithms that will be run in real-time by LIGO and 
VIRGO, enabling the observatories to distribute GW "inspiral alerts" to interested observers. 
The primary questions for the near future, then, are whether the GW inspiral signatures 
of these mergers will be detected, at what rate, and whether or not a coincident 7-ray or 
other electromagnetic signature can be identified. Unless our understanding is rather radically 
misplaced, all three questions are likely to be resolved once the era of the advanced GW 
facilities is fully under way. 

2.4- Choked Jet Supernovae 

The choked jet supernova, a theoretical scenario, represents the alternative fate that awaits 
a massive star if its core collapse generates a high-energy jet as in HL-GRBs, but in a fashion 
or within the context of a higher-mass stellar envelope that absorbs the jet energy before it 
is able to escape the star. As such, the chokcd-jct events evince no high-luminosity 7-ray 
emission; however, in the course of being quenched, their sub-stellar jets may undergo internal 
shocks which could accelerate protons and yield p7 and pp neutrinos in the TeV range [46-48] . 
Observationally, then, the choked-jet supernovae would be observed as a relatively low-redshift, 
highly-energetic supernova (a hypcrnova) with associated CN emission. 

Event rate estimates for choked-jet events are necessarily speculative. However, we note 
that if the process that leads to high-energy jet production during core collapse is largely 
independent of the nature of the overlying stellar envelope, then the choked-jet supernova rate 
may be comparable to the rate of HL-GRBs. Moreover, given that the coUapsar mechanism 
has multiple failure modes, including a high-mass stellar envelope, insufficient rotation, and 
too much jet precession [49], one can speculate that choked jet supernovae should, on generic 
grounds, be more common than HL-GRBs, although potentially (for these and related reasons) 
exhibiting less-luminous CN and/or GW emission. 



6 



2.5. Core Collapse Supernovae 

Core-collapse supernovae (CCSNe) are expected roughly every 30 to 100 years in our 
Galaxy; given the direct heritage of the revolutionary neutrino detection of SN 1987A [50], 
any such Galactic cataclysm would obviously represent a first-grade candidate multimessenger 
transient. 

Seen in relatively nearby galaxies (distances D ^ 100 Mpc), CCSNe are typically detected 
by optical observers within a few days, although as Swift has demonstrated [51], prompt de- 
tection of the high-energy shock breakout can be realized with sufficiently sensitive instrumen- 
tation. On this timescale, the prompt thermal {Ei, ^ 10 MeV) neutrinos should be detectable 
by IceCube and Super-Kamiokande for events within D ^ 50kpc [52]. 

CN may be produced by these "ordinary" CCSNe, either at shock breakout [41] or at 
later times, via interaction of the SN shockwave with a dense circumstellar medium [53]. The 
shock breakout CN emissions of ordinary CCSNe will be fainter than for the LL-GRBs, and 
so require very nearby events or larger-scale stacking analyses to detect with IceCube. CCSNe 
circumstellar interactions have the potential to generate a greater number of CN over an 
extended timescale of weeks to months, and may be detectable with IceCube, if candidate CN 
CCSNe can be identified with confidence and in sufficient numbers to make a sensitive search. 

GW emission from generic CCSNe are expected to be relatively weak [54], and hence, not 
detectable except in the case of a Galactic (or possibly, Small or Large Magellanic Cloud) 
event. 

Cosmic rays with energies below the "ankle" of the spectrum at E 4 x 10^^ eV are 
generally attributed to sources within the Galaxy [55]. Since no steady sources have been 
detected despite sensitive searches [56], they may be produced in occasional bursts within the 
Galaxy, including Galactic CCSNe. The decay length for relativistic neutrons is 9.2£^i8kpc, 
where E'lg is the neutron energy in EeV. Transient sources of neutrons of 1 EeV or above can 
thus be detected across much of the Galactic disk, including the central bulge region, and at 
higher energies also accessible with Auger, sources throughout most of the halo are detectable. 

2.6. Blazars 

Blazars are frequently detected in 7-rays at GeV (Fermi, AGILE) and TeV (HESS, VERI- 
TAS, MAGIC) energies; they are also highly luminous X-ray, optical, and radio sources. Since 
the non-thermal X-ray emission of blazars is generally understood to be synchrotron emission 
from electrons accelerated in shocks of an inner jet, shock-accelerated protons are expected 
as well. The target photons for interactions, leading to CN emission, could then be either 
synchrotron photons from co-accelerated electrons, or optical/UV photons originating in the 
accretion disk or scattered into the jet by broad-line region clouds. 

Blazars are highly variable, flaring sources that are much brighter across the EM spectrum 
during flare episodes than in quiescence. Expected CN fluxes for standard blazar models 
[21, 57, 58] suggest that typical individual flares cannot be detected with IceCube, with the 
summed contribution of ^100 such flares required to yield an expected >1 CN detection in 
IceCube. Separately, limits on the diffuse CN flux due to the summed contribution of all blazars 
in the Northern hemisphere have been published by IceCube using the 40-string (roughly one 
year integration) dataset, constraining some models [24], but without the statistical leverage 
that would be provided by comprehensive EM monitoring yielding the times and durations of 
blazar flares. 
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Individual spectral components of blazar flares are also variable, and the properties of the 
highest-energy components are poorly constrained due to the difficulty in obtaining sensitive 
TeV 7-ray data and simultaneous flare-triggered multiwavelength EM coverage. Moreover, 
recent Fermi and multiwavelength data suggest that spectral breaks and possible additional 
components are present during some blazar emission episodes [59-61]. Exceptional TeV 7-ray 
flares from blazars that lack associated X-ray emission, analogs of the lES 1959+650 "orphan" 
TeV flare [62-64], could be associated with hadron acceleration that would yield CN fluxes 
well in excess of those expected from typical flares, and hence, would be more readily detected 
by current facilities. 

Energetic GW emission from blazars in the frequency range of ground-based detectors is 
not anticipated. 

2. 7. Primordial black holes and other exotica 

If primordial black holes (PBHs) formed in the early universe with masses M ~ 5 x 10^^ g 
appropriate for them to Tuidcrgo explosive evaporation at the present epoch [65], then they 
would serve as a distinct and exotic type of multimessenger transient. 

As the black hole loses mass via Hawking radiation, dM/dt = —a{M)/G'^M'^ [66], its 
temperature Tbh increases, allowing an increasing number of particle types (degrees of freedom 
or dof) a{M) to be radiated. Once Tbh reaches the quantum chromodynamic energy scale of 
Aqcd ^ 200 MeV, free quarks and gluons will be emitted and fragment into hadrons, photons 
and leptons, resulting in a flux of high-energy cosmic ray neutrons and CN [67] that might be 
detectable from Auger and IceCube. 

Detection of PBHs would give dramatic confirmation of Hawking's hypothesis and theories 
positing early cosmological phase transitions (e.g., [68, 69]), and enable studies that would 
likely yield deep insights into ultrahigh-energy physics as well as quantum gravity. 

Other primordial relics, if they exist, might have decay modes yielding harder spectra than 
the PBHs (e.g., [70-76]). Power- law spectra with n oc E'"^ or harder would put out a majority 
of the total decay energy at the highest energies, meaning that the first detections, or most 
sensitive limits, on these processes are likely to arise from the highest-energy (largest area) 
facilities, IceCube and Pierre Auger. 

3. AMON Dataflow and Operations 

Multimessenger detection of one of the astrophysical sources described above, or of some 
entirely distinct cosmic phenomenon, will require coordination and cooperation between obser- 
vatories of strikingly different design and operation. These observatories produce data whose 
heterogeneity motivates the unified approach of AMON for detecting coincidences. Individ- 
ual facilities participating in AMON can be characterized as triggering facilities, follow-up 
facilities, or both. Triggering facilities are typically sensitive to one or more messenger type 
(photons, cosmic rays, neutrinos, or gravitational waves) and search in a wide field of view or 
monitor known sources for transient behavior. Prom their raw data, they construct candidate 
astrophysical events, here denoted "trigger events", which may transmitted through either 
public or private channels to AMON. These variously represent a single particle interaction 
{e.g., a muon neutrino interaction detected by IceCube), multiple detections processed into a 
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high level trigger {e.g., a Swift BAT light curve), or a sharp rise in a continuous measurement 
{e.g., an Advanced LIGO event). 

For AMON, a fundamental characterization of each observatory is its event- wise false pos- 
itive rate (FPR). The FPR is the dominant component of the total event rate that each 
observatory sends to AMON. The FPR is typically due to intrinsic detector noise or irre- 
ducible backgrounds that produce signal-like events. The FPRs for all triggering observatories 
considered here are listed in Table la. Later in this section we describe how these FPRs are 
calculated. 

By design, AMON will handle data streams dominated by these subthreshold signal-like 

events that individually cannot rise to the level of astrophysical discovery. For example, a single 
muon neutrino detected by IceCube could be of astrophysical origin, but the channel is indis- 
tinguishable from, and dominated by, atmospheric muon neutrinos. Similarly, single Galactic 
neutrons in Auger cannot be distinguished from the dominant background of charged particle 
cosmic ray events, and lowering the threshold of X-ray, 7-ray, or gravitational wave detectors 
increases signal acceptance at the cost of simultaneously admitting many additional accidental 
noise events. However, if two or more of these events are correlated in arrival direction and 
time they can attain a combined significance that results in AMON issuing an "Alert." An 
Alert either constitutes an immediate discovery or prompts follow-up observations that could 
demonstrate the presence of a "smoking gun" electromagnetic counterpart. By having all indi- 
vidual observatory FPRs on hand, algorithms running in the AMON framework are uniquely 
capable of calculating the aggregate Alert FPR associated with each coincidence detection. 
The Alert FPR is a statistical measure of the quality of a given multimessenger coincidence, 
endowing each Alert with a high degree of utility for multimessenger source searches. 

It is also possible for trigger events from a single observatory to be intrinsically above 
threshold, allowing the observatory to make a stand-alone claim for discovery. A burst of three 
neutrinos in close directional and temporal coincidence (IccCube or ANTARES), a burst of 
three air showers in coincidence (Auger), especially if their putative source is close enough for 
them to be neutrons, a statistically significant burst of hard X-rays (Swift BAT, Fermi GBM) or 
7-rays (Fermi LAT or HAWC), or a strong strain registered by a gravitational wave detector 
(LIGO- Virgo) would each constitute a discovery. For these stand-alone discoveries, AMON 
can serve as a quick and convenient conduit for disseminating the source coordinates, enabling 
timely follow-up observations at multiple disparate observatories. AMON will leverage the 
GRB Coordinates Network (GCN) [77], among others, to promulgate its Alerts. 

Underpinning AMON's sensitivity to new phenomena is the broad (multiple steradian) sky 
coverage and high duty cycles of the triggering observatories. Fig. 1 shows the results of a 
simulation of a calendar year for the triggering facilities that have signed the AMON MOU 
(IceCube, ANTARES, HAWC, and Swift BAT) and potential future signatories (Auger, Fermi 
LAT and GBM, and LIGO- Virgo), establishing that - absent severe disruptions - at least two 
facilities are always observing any given part of the sky, and that about 94% of 47r sr-yr lies 
within the field of view of three or more facilities. 

Existing multimessenger searches involving only pairs of observatories are typically in a 
master-slave relationship, e.g., a triggering observatory initiates follow-up at an optical tele- 
scope [78-80]. In contrast, AMON will enable archival and real-time searches for coincidences 
among multiple observatories in a peer-to-peer relationship that will provide markedly in- 
creased opportunity for discovery of new multimessenger sources, as detailed in Sec. 4 below. 



9 




Number of co-pointing observatories 




Right Ascension 



Figure 1: (left) Distribution of the total number of trigger facilities observing a cosmological source, 
averaged over time and sky location, (right) Average number of observatories simultaneously viewing 
a source, as a function of location. Spatial variance is caused by the limited field-of-view of some 
observatories and the movement of orbital telescopes. 



The real-time AMON analysis will be informed by the extensive effort already invested by 
the triggering observatories [78, 81-90]. The trigger events these observatories produce will 
be heterogeneous, covering the full range of messenger particles and differing significantly in 
background rate and localization precision. In some cases the events will include additional 
quality information that may be useful in refining the single-stream false positive rate. For 
purposes of AMON, all events will be couched in the common language of statistics, contain- 
ing information including a trigger time, event position, and a positional error or probability 
density function (PDF). 

Once the individual event streams have been characterized, including tracking observatory 
pointing and FPR, the multiple event streams can be combined. At the most basic level, 
this will be carried out via a coarse clustering analysis, searching for temporal and spatial 
coincidences. The output of the clustering analysis is dominated by pairwise coincidences, 
with the FPR for an ordered pair of triggers from observatories (a, b) given by^, 

4r-^^Ar(f]„,)Ai7,,, (1) 

where Ra is the subthreshold event rate from observatory a to AMON (assumed constant 
here), f^a is the field of view of observatory a, AT is the total width of the temporal search 
window, {Qab) is the time-averaged overlapping field of view of observatories a and 6, and 
AQab is the search area for coincidences between observatories a and b. In Table 1 (row "b"), 



^In some cases, it may be advantageous to use one event stream to trigger event selection for another, 
e.g. if the rate of single events is high, as with HAWC, or if the localization is poor, as with low-significance 
gravitational wave events. However, processing of an observatory's data can always be carried out, either at 
AMON or prior to transmission, to yield a manageable single-observatory FPR, and such preprocessing has 
been assumed for this discussion. 
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we show the result of this calculation for AT = 100 s and A^ab adjusted to give ~ 90% 
acceptance for signal event pairs (where each event is at the threshold for transmission to 
AMON). The self-coincidence of the individual observatories is given by the diagonal Raa 
(although this calculation is not appropriate in the case of preprocessed event streams) and 
the sum Rab+Rba {o, ^ b) gives the total number of pairs from a and b without respect to order. 
As expected, the combined FPR due to a pair of false positive events is orders of magnitude 
lower than that of the individual component subthreshold event streams. 

This particular calculation is by no means definitive. A trivial adjustment of the cuts will 
increase the rates (useful for generating AMON Alerts for testing purposes) . More importantly, 
additional information can drastically lower the FPR for the AMON analysis (Table 1, row 
"c"). For example, one can statistically stack the Alerts though an archival study, or shorten 
the temporal search window to search for exotic phenomena. In real-time, one can require 
that at least one event be of high significance or that the coincidence contains at least three 
events^. 

A promising feature of AMON is its ability to distribute Alerts in real-time to initiate 
follow-up observations. The relatively high FPR represented by the sum total of all pairwise 
coincidences among the observatories considered (Table 1, row "b") may be tolerable for some 
follow-up facilities, especially if their response is highly automated. However, many follow- 
up facilities will have limited observing time available for this science and may require a 
substantially lower FPR for triggering. This can be achieved via a refined analysis, providing 
a combined best fit location and error, as well as a maximum likelihood ratio or Bayesian 
probability. The likelihood or probability measure is drawn from a continuous distribution, 
in contrast to the FPR of the clustering analysis which, given a fixed event rate for each 
stream and a fixed set of cuts, takes on a unique value for any pair of event streams. Prior 
information can easily be incorporated into the likelihood analysis including, for example, 
previous detection limits, observatory sensitivity, or a galaxy catalog, and these choices can be 
tuned for different follow-up programs. 

As an independent channel of scientifically-relevant data and constraints, follow-up obser- 
vations have the potential to provide the added information required for a definitive discovery. 
By distributing candidate source positions in real-time to follow-up facilities, with statistically 
valid measures of the FPRs, AMON will enable efficient fast-response counterpart searches and 
studies across the electromagnetic spectrum. The existence and nature of such EM counter- 
parts may prove decisive in verifying the existence of some of the first multimessenger sources, 
and can reveal the nature and detailed properties of the source, including important ancillary 
information such as the source environment, host galaxy (if any), and properties of intervening 
gas and dust along the line of sight. The discovery potential of pairwise coincident events with 
positive follow-up detection is described in more detail below. 

AMON will be deployed in three phases. During the first year of operation beginning 
in mid-2013, AMON will use archived data to develop and tune analysis algorithms. Such 



^The FPR for a coincidence of the ordered triple (a, b, c) is given by iJ^^f ' ^ I §^ §t (ATA^abc)^ {^abc}, 
where the factors AQabc, an equivalent search area for a threefold coincidence, and (flabc), the joint field of 
view, are determined by Monte Carlo. The factor of | depends on how the temporal search is defined, where 
here we require that all three events occur within AT. In Table Ic, we present the totals li'^bc^ > 
case where all Auger data is included and for the case where Auger data is restricted to the galactic plane. 
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archival analyses, while a necessary first step prior to activation of real-time alerts, will also be 
of intrinsic scientific interest, and may point to possible new astrophysical signals or new and 
constraining upper limits on jointly-emitting source populations. Even as AMON moves toward 
real-time operation, certain source populations will be best discovered through a collaborative 
blind study of accumulated data; for example, in searches for exotic particles, as discussed 
below for primordial black holes. To enable this, participating AMON members will, at any 
time, be able to utilize the AMON database and data products, under rules that are established 
in the AMON MOU [91]. 

In subsequent operational phases, standard AMON analyses will be modified to operate 
in real-time, enabling identification of candidate transient sources as soon as the events are 
transmitted. At first, the effort will focus on those discoveries that can be made by the 
triggering facilities alone, with FPRs similar to those in Table Ic. Even at this stage, the 
real-time nature of the network plays a critical role, providing the triggering facilities with 
AMON Alerts that can prompt a closer analysis of their own data. The Alerts shared at this 
stage will include many whose combined FPR is too high to claim detection of a transient. 
These will not only exercise the AMON system, but may initiate a deeper search by those 
participating observatories that did not initially report a triggering event. An example of 
this discovery mode is the joint search for high-energy cosmogenic neutrinos and gravitational 
waves, as discussed below. 

As quickly as is feasible, AMON will transition to its final phase, with Alerts distributed 
to participating follow-up observatories to enable the near real-time search for electromagnetic 
counterparts. For sources that exhibit detectable afterglow, this is the most powerful technique, 
leveraging the high acceptance of the initial FPR threshold, yet resulting in a final FPR that 
is lower than that available by the other techniques. The estimated final FPRs are shown in 
Table Ic, for the case of serendipitous follow-up detection of GRB afterglow (assuming Swift 
XRT sensitivity and a ~ 1° search region) or a SNc light curve (assuming the detection limits 
of [78]). In the next section, as a proof of concept, we show that significant sensitivity gains 
are made by AMON in the search for electromagnetic counterparts to cosmogenic neutrino 
sources. 
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4. Discovery Power of AMON 

With the array of sources and messenger types described in Sec. 2 and the possible groupings 
of observatories described in Sec. 3, there are a multitude of ways in which AMON may provide 
scientific gain. Unable to cover all features of AMON in one paper, we focus here on three 
examples. Sec. 4.1 examines how the distribution of AMON Alerts can greatly enhance the 
search for electromagnetic counterparts to candidate sources of cosmogenic neutrinos, while 
Sec. 4.2 quantifies the gains for a search for joint sources of gravitational waves and CN. 
Sec. 4.3 discusses an example search for exotic phenomena, showing how an AMON analysis 
may provide the first clear signature for primordial black holes. 

4.1- Follow-up of candidate cosmogenic neutrino sources 

As a first example, we consider the follow-up imaging of candidate cosmogenic neutrino 
sources. Such searches are already underway via efforts by individual neutrino observatories 
and their follow-up partners. For example, in [92], the IceCube Observatory identifies clusters 
of two or more neutrino events as a trigger for follow-up imaging. Due to the high rate of 
atmospheric neutrinos and other backgrounds, the expected FPR for a pair of up-going events 
for IceCube-86 is approximately 30 yr^^ (Table 1, row "b"). Alerts are then distributed to the 
ROTSE telescope array or, after using a maximum likelihood analysis to refine their number, 
to Swift and PTF [80, 93]. This program is already producing the first limits on SNe with 
choked jets [78]. 

As a shorthand, we introduce the notation "2z/" to indicate a trigger condition of 2 or more 
coincident neutrinos, or v-Nj for one or more neutrinos and N or more 7-rays. To be concrete, 
we study the follow-up imaging of 2i^ alerts at a rate of 10 yr~^ using Swift's co-aligned X-ray 
and UV/Optical Telescopes (XRT, UVOT), with this considered a status quo approach for the 
scientific community. Due to the narrow field of view (0.4° diameter) of the XRT and UVOT, 
7 Swift pointings are required to cover the error region of each alert, totaling 70 pointings in 
all. 

In contrast, the AMON approach enables the realtime coincident analysis of neutrino and 
electromagnetic data streams, comparing the arrival time and direction of a single candidate 
CN with data from Swift's large field of view Burst Alert Telescope (BAT), Fermi's Large Area 
Telescope (LAT), or the upcoming ground-based HAWC. Because of the superior AMON lo- 
calizations, derived via the joint analysis with high-energy 7-rays, the same follow-up telescope 
resource (70 pointings) can be used for follow-up of a significantly-increased number of AMON 
alerts. Here we consider a program that follows up five 2u alerts (at 7 pointings each) and 35 
u-N^ alerts (1 pointing each). For the u-Nj alerts, we considered the pairing of IceCube with 
both BAT and LAT. 

The relative sensitivity of the status quo and AMON approaches was studied via a Monte 
Carlo calculation. As a source model, we assumed that both 7-ray and neutrino spectra follow 
a broken power-law with initial and final logarithmic slopes of a = 1.0 and /3 = 2.0, and break 
energies £^ = 0.2 MeV and = 0.35 PeV, respectively. The neutrino and 7-ray fluence was 
allowed to vary, as indicated by the axes of Fig. 2. Using realistic models for the IceCube-86, 
BAT, and LAT effective areas, point spread functions, and their overlapping field of views, we 
were able to determine the acceptance of each simulated source. Realistic background models 
were adopted for each observatory, enabling the calculation of false positive rates for each joint 
analysis. 
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In each case, the FPR was tuned via a hkehhood analysis of the relative position of the 
V and 7 events. Subthreshold BAT and LAT events were allowed into the simulated data 
stream, with terms added to the likelihood function to favor the more significant EM signals. 
Furthermore, we added a regulating term to the likelihood function to account for the possi- 
bility of a real 7-ray source in accidental coincidence with an IceCube background event (e.5. 
an atmospheric neutrino). Once the FPR is chosen for any pair of observatories, we identify 
the required likelihood threshold, and can then calculate the corresponding sensitivity to our 
modeled source. While the initial FPR is quite high (10-70 alerts yr~^), the final FPR will be 
orders of magnitude lower after the inclusion of follow-up data and, as such, can be ignored 
in the study of sensitivity. While here we have characterized the sensitivity of the triggering 
methods (AMON versus status quo), there are a number of unknown factors in the sensitivity 
of the subsequent follow-up searches; for example, due to optical magnitude. However, we 
study the ratio of the sensitivity of the two methods, expecting the acceptance of the two 
follow-up programs will be approximately the same and to cancel in the ratio. 



Figure 2: Sensitivity of the AMON multi-channel 
method compared to the status quo, presented as 
a ratio ranging up to ^ 10^, for follow-up observa- 
tions of candidate astrophysical neutrino sources. 
The assumed shapes for both gamma and neu- 
trino spectra follow a broken power law with ini- 
tial and final logarithmic slopes of a = 1.0 and 
/3 = 2.0, and break energies £.y — 0.2 MeV and 
£1, = 0.35 PeV. The vertical line shows a typi- 
cal threshold for a hard X-ray observatory {e.g., 
Swift BAT), and diagonal lines corresponding to 
predicted neutrino to gamma-ray fluxes under a 
range of current theoretical GRB models [22] . 

The results are shown in Fig. 2, where substantial gains in sensitivity (up to 10^) are 
realized. The sensitivity of each method is defined as the fraction of true sources pursued 
using a given follow-up capability. In addition to enabling pursuit of an increased number 
of alerts overall, the AMON approach realizes these large sensitivity gains by triggering from 
single neutrinos rather than a neutrino pair, leveraging the increased interaction strength of 
7's over j/'s. The left axis of Fig. 2 is defined in such a way that it could equally well describe 
the expected event rate for ANTARES, with only a small change in the contours due to the 
narrower neutrino point spread function. 

Gains of this type have also been demonstrated for gravitational wave source searches, 
using GW-7 triggering over the status quo GW-only methods [94]. For GW-only searches, the 
poor GW localization is filtered by requiring coincidence between the GW localization and a 
nearby cataloged galaxy [79, 95]. While useful, this approach suffers from several shortcom- 
ings: First, it depends on the completeness and astrophysical relevance of the associated galaxy 
catalog; second, it cannot be extended to the subthreshold regime due to the dramatic increase 
in false alarm rates; and finally, the approach is limited to an horizon of ^100 Mpc, beyond 
which galaxy fields become crowded and catalogs are highly incomplete. The multi-channel. 




10^^ 10"^ 10° 10' 10^ 
X-ray fluence in 15-150 keV [N.cm"^] 
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multimessenger approach by contrast provides a unique candidate source position for follow- 
up pointing, and leverages the coincidence technique to push into the subthreshold, possibly 
revealing otherwise-undetectable source populations. With Advanced LIGO expected to ulti- 
mately achieve sensitivity out to ;^500Mpc [96], the AMON approach can offer a >125-fold 
gain in sensitivity for events with bright EM counterparts, or > 1000-fold for regions that have 
the least complete galaxy catalog. 

4-2. Joint detection of neutrinos and gravitational waves 

Here we consider the joint detection of gravitational waves and high energy cosmogenic neu- 
trinos as an example of a science investigation enabled by the realtime operation of AMON. 
The search for jointly emitting GW+CN sources has been addressed in [97, 98], including pre- 
dictions for Advanced LIGO- Virgo and IceCube-86, and we follow the same general approach 
(see also [99, 100]). However, while previous authors focus on the joint analysis of a multi- 
messenger data archive, we emphasize the possibility of carrying out this analysis in realtime. 
We characterize the sensitivity of a blind (untriggered) search by Advanced LIGO- Virgo using 
the horizon distance -DgW; inside which a GW source is expect to generate an above-threshold 
detection. 



where is the isotropic equivalent gravitational wave energy of the source. This is based on 
the results of [101] and assumes little of the GW waveforms. Without AMON-likc infrastruc- 
ture, a real-time detection of both GW and CN requires that a source be within this horizon 
and that it produce a neutrino flux sufficient to trigger > 3 events in IceCube or ANTARES 
(since rij^ < 2 would lead to an inconclusive false positive rate) . These detections would then be 
compared post facto, either through an archival analysis or above-threshold alert distribution 
{e.g. through GCN), to determine if there are any joint GW-CN sources*^. 

By contrast, a system like AMON will provide two-way communication between observa- 
tories, enabling a realtime triggered search of GW data, either by transmission of pre-alerts 
from AMON to LIGO- Virgo or by carrying out a triggered search on subthreshold GW data 
transmitted to AMON (albeit in some preprocessed form) . Searches of this kind suffer a lower 
trials penalty and may assume various templates for the sources. For example, the postulate 
that CN are detected in coincidence allows one to make an assumption about the orientation of 
the source, making the GW signal more favorable for detection [102]. In all, we estimate that 
these effects will combine to increase the horizon distance to 2 x Dqw [102]. Furthermore, only 
a single neutrino event is required to trigger the analysis, providing significant gain for this 
approach since the expected number of neutrinos has been observed to be typically ^ 1 [103]. 

Since we are considering the relatively deep horizon distance for Advance LIGO- Virgo (so 
that galaxy distribution can be assumed uniform), but not so far that redshift effects become 
important, it is sensible to apply the effective volume approach of Appendix A, integrating up 
to the physical limit of Dqw or 2Z)gw- The resulting effective volumes are given in Table 2. 



''in [97], the authors consider a shghtly different problem, being interested in the above-threshold detection 
of GW or CN. While an interesting measure, their calculations do not provide the sensitivity for above-threshold 
detection of joint sources that we need for comparison to AMON. 




(2) 
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Table 2: Effective volumes for a joint GW-CN search, 
assuming E'^'^ = IQ-'^Mq and E'J° = 10"*^ erg for a 
choked GRB model. 



Analysis 


Integration 


Neutrino 






limit (Mpc) 


threshold 


[Mpc^] 


GW+CN 


80 


ri/y > 3 


4.6 X 10^ 


AMON 


160 


ny>l 


5.9 X 10^ 



Under the assumption of a low false positive rate, the 90% upper limit (chosen for consis- 
tency with [97]) for the joint GW-CN source density is given by 



Pgw-CN < 



2.3/b 



[yr~^Mpc 



-3i 



(3) 



where T is the livetime, which we take to be one year, V^^^^ is the solid angle averaged effective 
volume for the given method, and is the beaming factor for neutrinos. We take the latter 
to be fb = 14 (from [97]) but it will cancel in our final result. 

The effective volume depends on -Eq^ through equation 2. For neutrinos, the expected 
number of detections from a fiducial source at distance r can be related to the neutrino isotropic 
energy -BJf", using a linear relationship provided in [97], 



El 



if- 

K VlO^Serg 



10 Mpc^ 



(4) 



where, for IceCube-86, k = 0.75 for a high luminosity GRB model and k = 1.6 for a model of 
choked GRBs. 
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Figure 3: (Left) Projected GW-I-CN source population upper limits after one-year of AMON utilizing 
observations with Advanced LIGO- Virgo and IceCube-86. (Right) Ratio of the sensitivity for the 
AMON approach to the sensitivity of the status quo approach (as described in the text). 
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The results are shown in Fig. 3. The left side shows the upper limit in units of yr~ Mpc~ 
that would be obtained by a non-detection with AMON. The right side shows the multiplicative 
gain in sensitivity that the AMON approach provides over an analysis that requires above- 
threshold detection in both GW and CN channels. The gain is most significant for bright GW 
sources and weak neutrino sources, where the burden for discovery is shifted from the CN to 
GW channel. However, the gain is strictly greater than 1, with an asymptote of 8, representing 
the increase in effective volume due to triggering the GW search. 

4-3. Searching for primordial black holes 

Here we consider the search for primordial black holes as an example of a collaborative study 
of accumulated multi-facility data, as enabled by the AMON infrastructure and data sharing 
policies. In particular, we have modeled a joint search with the IceCube neutrino observatory, 
HAWC TeV gamma-ray observatory, and the Pierre Auger cosmic ray observatory. All are 
sensitive to particles produced in the final stages of PBH evaporation, including a potential 
signal from ultrahigh-energy neutrons detected at Pierre Auger, since these neutrons would 
not suffer from the magnetic deflection and time delay effects of charged cosmic rays (the chief 
UHE neutron background). 

Based on the energy threshold of these experiments and energy dependence of various 
models predicting energy spectra of PBHs in their final stage of evaporations {e.g. [104-106]), 
the model of a non-rotating, uncharged black hole without a chromosphere by MacGibbon 
and Webber [104] was chosen for estimation of the expected signal at each of these detectors. 
The time integrated particle spectra above 100 GeV were computed following methods from 
refs. [107, 104, 108, 109]. The main particle decay chains considered in our calculations are 
described in Appendix B. 

It is anticipated that, when fully operational, HAWC by itself will either detect or provide 
the best upper limit on direct searches for PBH decay, around 2 orders of magnitude better 
than the current upper limit (as shown in Fig. 4). To do so, HAWC must perform a blind 



search of its own data. After a cut to remove many hadronic events, 5.1 kHz of subthreshold 
trigger events remain, which will be analyzed for coincidences across a spatial template with 
10^ trials [110]. Taking a generous time window of AT = 5 s and a number threshold of n > 20, 




10^ 



10^ 



10' 10' 
PBH explosions (pc"^yr') 



Figure 4: Expected number of sources detected in at least two 
of the v, 7, and n channels, modeled here with the IceCube- 
86, projected HAWC, and Auger sensitivity, as a function of 
the primordial black hole density. The shaded region to the 
right indicates PBH densities that previous direct searches 
have excluded, while the vertical line shows where the up- 
per limit is expected to move if HAWC does not observe a 
gamma ray excess after a year of livetime. In between, there 
is the possibility of positive joint detections, with potentially 
a unique PBH signature. 
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the HAWC false coincidence rate can be reduced to < 1 yr~ . However, it is expected that 
HAWC will also observe real astrophysical transients that are unrelated to PBHs, forming a 
second type of false positive for the PBH search. We conservatively estimate this number to 
be -Rhawc < 17 bursts per year, approximately equal to the known LAT GRB rate, allowed to 
fluctuate up by 2cr. If, after T = 1 year of livetime, 17 unidentified bursts were observed, the 
99.73% (3(7) upper limit for PBH would be given by 

32 1 

PPBH < [yr-^pc-3] (5) 

where 32.1 is chosen to give the 99.73% upper limit of the Poisson distribution, and V^^^ is the 
integral of the trigger condition P{n > 20\r) over Anr'^dr, averaged over solid angle, where r is 
the distance from the Earth. The value of V^^^"^ is specific to the sensitivity of the observatory 
and the intrinsic 7-ray fluence of the fiducial source, and is further described in Appendix 
A. For the current calculation, we estimate V^^^ = 0.037 pc~^. The calculated value of 
Ppbh < 3.5 X 10^ yr~^pc~^ can be compared to the current best limit of 1.9 x 10^ yr^^pc"^ 
from VERITAS [111]. 

Between the current VERITAS limit and the projected HAWC limit, there exists the pos- 
sibility of a future positive detection of PBH decays. However, the HAWC data by itself may 
be unable to distinguish between PBH and other scenarios. A joint AMON study of HAWC 
data with neutrino and neutron data (modeled here with IccCube and Auger), will enable 
the search for multimessenger coincidences, providing a unique PBH signature. Just a few 
such coincidences should reveal a timing and energy structure that are indicative of Hawking 
evaporation. 



Table 3: Effective volumes for multimessenger PBH search. Values are 
averaged over Att sr, with only those regions of overlapping sensitivity 

contri1)uti]ig lo each i)air\\is<' calculation. 



Observatories 


Tri 


gger 


conditions 


T/(cfi) 




^SUSY 


A 


B 




A 


B 


[pc3] 




[pc3] 


HAWC 




ri'y 


> 20 




0.0374 




0.245 


IceCube 


HAWC 




> 1 


ra-y > 13 


9.8x10- 


■5 


6.4x10-^ 


IceCube 


HAWC 


rii, 


> 2 


n-y > 13 


1.3x10" 


-6 


8.3x10-6 


HAWC 


Auger 




> 13 


n„ > 1 


1.3x10" 


'7 


8.6x10"^ 


Auger 


IceCube 


rin 


> 1 


n„>l 


4.9x10" 


-9 


3.2x10"^ 



Folding together the PBH source model from Appendix B and the volume integrals of 

Appendix A, we calculate the effective volumes given in Table 3. To maximize sensitivity 
to neutrinos and neutrons, we use number thresholds of rij^ > 1 and n„ > 1. For 7-rays, the 
low false alarm rate can be relaxed significantly by requiring coincidence with a single v or n, 
allowing the HAWC threshold to be lowered to n-y > 13. The expected rate of true positive 
joint detections in channels a and h is then given by 

= [y^-Vc-] (6) 
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which we show as a function of ppBH in Fig- 4. Over the range of p which might lead to future 
positive detections by HAWC, the total number of expected multimessenger coincidences ranges 
from 0.1 to 100 yr~^. Most of these are v-N^ {i.e. a coincidence between a single neutrino 
and N 7-rays). However, there is a non-negligible possibility for observing n-N'y coincidences 
(where the neutron is detected by Auger) or 2iy-N'y, providing an exciting possibility for 
detecting three subthreshold events in coincidence (where the Nj from HAWC are treated as 
a single subthreshold event). 

5. Summary and Conclusions 

The Astrophysical Multimessenger Observatory Network under development at Penn State 
will link multiple high-energy, multimessenger, and follow-up observatories together into a 
single comprehensive system. We have explored the scientific opportunity for AMON (§2), 
which is centered on the discovery and exploitation of multimessenger transients. The nature 
of the brightest such transients, which may manifest as the first detections of gravitational 
waves, high-energy cosmogenic neutrinos, or high-energy cosmic ray neutrons, is still unknown; 
likely candidates include blazar flares and a diverse array of 7-ray bursts and supernovae, as 
well as exotic phenomena such as the evaporation of primordial black holes from the early 
universe. 

We have described the design, infrastructure, and current and projected partner facilities 
of AMON (§3), showing that the wide fields of view, high duty cycles, and subthreshold event 
rates of the facilities are such that a robust and automated statistical search for coincident 
events seen in the data of two or more facilities is both interesting and feasible (Table 1). 
In particular, false positive rates for the resulting AMON Alerts are sufficiently low that 
comprehensive ground-based optical follow-up campaigns can be contemplated. 

With neutrino and 7-ray observatories are already participating in AMON, and discussions 
underway with cosmic ray and gravitational wave observatories, as well as multiple follow-up 
observatories, AMON is poised to begin real-time operations within a year. 

To demonstrate the power of the AMON approach, we have carried out three sets of 
simulations using current theoretical models of multimessenger phenomena (§4). In the first 
simulation, we explored the challenge of EM counterpart searches for candidate GRB-associated 
CN observed in IceCube, showing that vetting candidate CN events against multiple EM data 
streams realizes a > 1000-fold gain in the efficiency of EM follow-up observations (Fig. 2). In 
the second simulation, we explored the improvement in search sensitivity for jointly emitting 
GW-I-CN transients that is realized by extending these searches into the subthreshold regime, 
rather than restricting the search to events that generate statistically-significant signals in both 
channels; we find that a > 10- fold increase in event rates, or improvement in upper limits, is 
easily achieved (Fig. 3). In the third simulation, we explored the multimessenger signature of 
PBH evaporation, demonstrating that - consistent with current upper limits on the local PBH 
density - the coincidence of an IceCube-detected neutrino with a cluster of HAWC-detected 
7-rays could provide evidence for observation of a PBH evaporation event within the first year 
of HAWC operations (Fig. 4). 

In addition to carrying out the real-time multi-facility transient searches needed to realize 
these and other science gains, and enabling follow-up EM observations by distributing transient 
alerts to interested observers, AMON will provide a framework for large observatory coUab- 
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orations to work together on targeted archival analyses using AMON's comprehensive events 
database and coincidence analysis toolkit. In these several ways, AMON will leverage and 
enhance the capabilities of existing and future high-energy and multimcsscnger observatories, 
powering a new and ambitious exploration of the transient universe using all four forces, and 
so helping to realize the immense promise of this dawning age of multimessenger astronomy. 
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Appendix A. Effective Volume Integrals 



In a number of examples in the text, a simple population model was applied that assumes 
a uniform spatial distribution of sources with a fiducial set of source parameters (encapsulated 
here by the parameter A). In the case where z ^ 1, we can define an effective volume 
for a specific set of trigger conditions by the Euclidean volume integral, 



[ 47rr^drP(true|A,r), (A.l) 
Jo 



where P(true|r, A) is the conditional probability of a true positive detection, given the fiducial 
source parameters and distance to the source r. The maximum distance D may represent 
either a physical limitation of the method or else a parameter that we use to regulate the 
integral before taking D — ^ oo. The parameter (s) A may depend on pointing direction, and so 
we have written V^'^^^ as difi^crcntial with respect to solid angle ^l. The expected rate of source 
detections is then given, after integrating over solid angle, by 

ii; = pF(eff), (A.2) 

where p is the intrinsic rate per unit volume of transient sources. In the special case where the 
observatory is equally sensitive across its field of view, ^^^^p- is independent of pointing and 
we can write 

R = pn—. (A.3) 

In the case where the participating observatories are each monitoring for Poisson processes, 
P can be taken to be the cumulative Poisson probability distribution above some set of number 
thresholds no- The vector notation includes the possibility of studying events from multiple 
observatories, with indices a = 1 — >■ M, so that 

P(n>no|A,.) = n ' 

a=l \ ka=0 «• \ / y 

where Xa = ^ J dE^^Ai^^^ (E) is a measure of the number of particles that can be detected by 

observatory a (where Na is the total number of particles produced and a'^^ is the observatory's 
effective area). It is normalized in such a way that Xa/r'^ is the expected number of detections 
for a source at a distance r. The trigger condition n > ng is shorthand for simultaneous 
detections above the specified number threshold for each observatory. 
As such, we are interested in integrals of the form 




-^A:(A)r{K-^,Aj, (A.5) 



22 



where we have changed the integration variable to x = A/r^ and defined A = Yla^a, = 
^„ kai and the term Qg(A) = ^^^-^ — ^ Ha (x)*^"- multinomial theorem, we note that 

where we have used a single subscript on the left hand side (rather than a vector). As expected, 
this is the same result as equation A.5 applied to a single observatory. The interpretation is that 

equation A. 6 can be applied to triggers from multiple observatories that are interchangeable 
{e.g. when considering k = 2 neutrinos form IceCube or ANTARES, regardless of which of the 
observatories triggered) . 

Importantly, vq and vi diverge as D ^ oo, but the linear combination ^ — vq — vi is finite. 
Thus, if we are looking for a single Poisson trigger, then a physical limit D must be applied to 
achieve a finite result {e.g. by requiring that a single detected neutrino originate from within 
the sensing region of a gravitational wave network). However, if we instead require that n >2, 
then the result is finite as D — >■ oo. 

We can also marginalize over one observatory, defining k = {k\,k') and using the normal- 
ization of the Poisson series to find 



ki=Q 

We can apply this result to two observatories. In the case where we search for at least m events 
from the first observatory and at least n events from the second, 

jT/(eff) r)3 

In- = ^ -(M)(Ai,A2) 

i,j>m,n 

^3 m— 1 n—1 m— In— 1 

= ^-E^^(^i)-E^^(^2)+ j;^;(,,i)(Ai,A2). (A.8) 

i=0 j=Q i=0 j=0 

However, if we search for a total of n events, regardless of which observatory they come from, 

j=0 

Both equations A.8 and A. 9 lead to finite results as D — > oo. 

We may wish to consider the case where one or both observatories have an effective area 
that changes significantly off-axis, meaning that ' depends on direction. In general, one 
must carry out the integration over solid angle numerically, although we have already noted the 
trivial analytic result in the case of equation A. 3. In addition, there is an analytic solution in 
the case where only the first observatory has a dependence on zenith angle 0z , that dependence 
is proportional to cos^^, and the number threshold of the first observatory is ni > 1. However, 
the result is sufficiently complicated that we choose to omit it here. 



JX,,X') = vp{X'). (A.7) 
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Appendix B. Primordial Black Hole Source Model 



Here we describe the main decay chains leading to the final particle spectra from the PBH 
explosion under the assumptions of the Standard Model (SM) and Supersymmetry (SUSY). 

Appendix B.l. SM decay chains 

During the final stage of the PBH life, all of the 118 SM particle degrees of freedom {dofj 
are radiated away. Since there are more dof for quark and gluons (72) than for leptons and 
photons (26), the final spectra of both HE and UHE gamma rays and neutrinos originate 
mostly from the decaying hadrons [104]. Neutrons and antineutrons are produced from these 
decaying quark fragmentation products, as well. 

We assumed that quark and gluon jets emitted during the final stage of PBH evaporation 
fragment into pions (vr^, 7r~ and vr"'') and baryons (p, p, n, and n) with branching ratios of 0.97 
and 0.03 [112], respectively. A fragmentation function of quarks or gluons into a particle 
of type X has a form [112]: 

%-^-\l-x)\ (B.l) 



Ax 16 

where x = Ex/E, Ex is the energy of a particle of type X, and E is the quark/gluon jet 
energy. 

For baryon spectra, we assumed that equal numbers of p, p, n, and n are produced after 
hadronization. Thus, the final n+n spectrum is obtained by multiplying distribution by 
0.015, convolving it with the Hawking primary spectra for quarks and gluons and integrating 
over time. 

The photon spectrum is obtained by convolving a flat photon energy distribution from 7r° 
decay with the 2-1/3-0.97- (where factor 2 accounts for the number of photons produced 
from each decaying pion and 1/3 is a fraction of neutral pions to the total number of pions) 
and the Hawking primary spectra. 

Prom decaying tt" and tt"^, and Vf^ are produced, together with /x~ and /x^. These 
neutrinos have fiat energy distribution, but there are also secondary neutrinos originating 
from the consecutive muon decays: iJ,~^^~^ — t- e"''^"-' + i^eii'e) + J^(t'^)- Energy distributions of 
each of these neutrinos are convolved with 2/3-0. 97- and the quark/gluon Hawking spectra, 
and summed afterwards. Since production of from decaying hadrons is highly suppressed, 
most of this neutrino flavour originate from the direct Hawking radiation and from the decays 
of the directly emitted r leptons. Therefore, one would expect a flux of i/^ that is almost two 
orders of magnitude lower than the total flux of other neutrino flavors. We neglected these 
two Vr contributions to the total neutrino flux in this work. 

Appendix B.2. SUSY decay chains 

If the Minimal Supersymmetric Standard Model (MSSM) describes Nature at the high- 
energy scales associated with the explosions of PBHs, the number of do/ available to be radiated 
increases by more than a factor of two: 244 compared to 118 dof in the SM, where we neglected 
graviton and gravitino c?o/ and included five physical Higgs fields (see, for example, Rcf. [113]). 
This leads to an increase of the factor a(M) that directly affects the rate of evaporation by a 
factor of ^ 3.3 when the PBH temperature reaches the SUSY particle production energy scale. 
Prom that point, the time left until complete evaporation will be shorter compared to the time 
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predicted under the SM assumption. For example, one second before complete evaporation, 
the temperature of the MSSM PBH would be T ~ 5.5 TeV with almost 50% more mass to 
be radiated than in the case of the SM PBH with T ^ 8TeV at the same time left before its 
complete evaporation. 

Given a number of unknown supersymmetric parameters (more than 100), for simplicity, 
we assumed that gluinos (g) are heavier than squarks (q) (as in, for example, the SPSla 
benchmark scenario [114] and the mSUGRA B benchmark model [115]), thus they decay into 
antisquark/quark and squark/antiquark pairs [116]: 

g^q + q,q + q. (B.2) 

The right chiral states of squarks would decay then mostly into the lightest netralino (x?): 

^ X? + q, (B.3) 

whereas the left handed states would decay into charginos (xf ) or heavier neutralinos (X2) [117]: 

qL^xf + Q, X2 + 9- (B.4) 

Further down these decay chains, we assumed that xf and X2 preferably decay into leptons 
and charged leptons plus the lightest neutralino, respectively: 

xt ^ i^^, X2 ^ i^rxi- (B.5) 

Without carrying a complete calculation, our rough estimate for the multimessenger particle 
fluxes under the assumption of the MSSM are: ^^4 of the SM flux for neutrinos and '-^3.5 of 
the SM flux for photons and neutrons/antineutrons. 
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